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Abstract 

Coherent dynamics of two-, three-, and four-leVel quantum 
systems, simultaneously driven by concurrent laser pulses of 
arbitrary and different forms, is treated by using a 

nonperturbative . group-theoretical approach. The respective 
evolution matrices are calculated in an explicit form. General 
aspects of controllability of few-level atoms by using laser 
fields are treated analytical ly. 

1. INTRODUCTION 

We analyze the general aspects of the problem of dynamical 
coherent control of atomic populations in the framework of 
semiclassical approach when an external C laser! field is 
considered as an inexaustible energy reservoir and only the atomic 
behaviour is considered to be controllable. The natural 

controllers are strengths, frequencies and phases of the 
components of a polychromatic laser field driving simultaneously 
different atomic transitions. An initial atomic state may be also 
considered as a controller. 

The dynamics of any nonstationary quantum system is described 
by the time evolution equation 

i ^ UCt,01 * HCt) UCt,0) , UC0,0! * I , h ■ 1 Cl! 

with the Hamiltonian 

m 

HCtl = Z h. Ctl H. . C2I 

k=1 * * 

which depends on time explicitely. It is not an easy task to solve 
the eq. 1 analitically even for a two-level system exposed to an 
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arbitrary varying external field. Usualy used adiabatic and 
weak-field approximations as well as the approximation with slowly 
varying amplitudes breaks down if we deal with short and/or 
intense laser pulses. The well-known Floke method is useful only 
in the case of periodic excitations. The high level of control of 
atomic populations by using designed laser fields calls for a new 
nonperturbative 'technique that should be, on the one hand, rather 
general to enable one to treat uniformely the atomic dynamics 
under various excitation conditions, and on the other hand, it 
should be able to result in explicit solutions of the evolution 

equation C1T. 

The group-theoretical technique will be applied to solve this 
task for two-, three-, and four-level atoms exposed to an 
arbitrary time-varying polychromatic laser field. It is based on 
the concept of dynamical symmetry of nonstationary quantum 

processes fl-31. 

After solving the eq. 1 for a given model system, the problem of 
finite control 141 is to find such values of controllers Th^CtD, 
j*l,...,m> that enables us to transform the system from an initial 
state to the desired final state at the target time T > 0. Such a 
transition is governed by the equation 

| sCT) > = UCT, 01 | sCOI >. (31 

Generally speaking, the solution of this task is not unique. In 
addition one can put the task of optimal control which is to find 
such values of controllers that minimizes a given loss functional 
or maximizes a desired quality criterion [51. 

2. GENERAL BACKGROUND 

A single atom with N nondegenerate and nonequidistant levels 

H o IV = w » IV> m ' C45 

interacts with amplitude- modulated components of a polychromatic 
f ield 

N N 

E = Cl/21 J 1 E kl fU exp Ci w kl U + c,c * C5!) 

1> k 
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by such a way that only one fk, 1) -component of the fifcld Ej,j is in 
resonance with the Ck, 13- atomic transition. Tf one deals with the 
electric dipole transitions 

H int * - FXU rtru • 

then by the parity reasons the energetic matrix Hq + H in ^ turns 
out to be tridiagonal. However, we will treat the general case 
Cwhen all the NCN-D/2 atomic transitions are driven independently 
on each other!) to be able to include into consideration, if 
necessary, other interactions that may couple nonadjacent levels. 
The complex variable E^Ct) measures amplitude, phase, and 
polarization of the respective component. 

Writing the state vector of a N- level atom in the form 

ism> = L c « cu 'v exp r_i w » *•’> f65 

m=l 

one can obtain from the nonstationary Scrodinger equation the 
following set of equations for the probability amplitudes in the 
rotating wave approximation 

N 

*V tJ ■ y h ki ft/) s i fu + c * r *» cn 

where 

h kl CtJ s " n/?J E ki m d ki m * f 8") 

The resonant Hamiltonian has the following structure 

H nn m = °> H kl CU = h kl ftl 88 H Tk Ct) » k * u 

Finally, the probability of finding our N-level system on level 
n at time t is expressed in terms of the evolution-matrix elements 
and the initial populations as follows 

P n c» ■ l , l%tn | 2 p.coj. C10J 
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3. COHERENT CONTROL ON THE SUC21 DYNAMICAL GROUP 


A quantum system with the SUC21 dynamical symmetry is a 
fundamental model in the semiclassical theory of field-matter 
interactions. The herraitean Hamiltonian of such an arbitrarily 
driven system can be cast in the form 

HCO * h Q Ct1 Rq + h*Ct1 R_ + hCtl R + C11I 

with the generators satisfying the commutation relations 

IR + , R_ 1 * 2Rq, l R q , R +> J * ±R +> _ , C12:) 


where c-number parameters h Q and h are assumed to be arbitrary 
analytic functions. Writing the evolution operator in the 
factorized Wei -Norman form 161 

U * exp Cg 0 - iJg h Q CTl dTl Rq exp g_R_ exp g + R +> C131 


one can obtain the governing equation for g * expCgg/PJ 

g’ - C £ + i h Q ) g + |h| p g * 0, gCOI * 1, gCOI * 0. CHI 

Once the eq. 14 is solved all the other SUC21 group parameters are 
found in quadratures 16,2,31 It follows from the unitariry of C131 
the simple conservation laws 

|g| 2 - g_ g + s Igl 2 c l + |g + | 2 i * l. cisi 

In the two-dimensional representation we get the expression 


Cl/21 

U 


g 12iCg1* + h 0 g*l/2h* ' 

C2i g -hgg1/2h g * 


C16I 


which is the time-evolution matrix for a two-level system driven 
by an external field with arbitrary amplitude and frequency 
modulation. In terms of populations the complete formal solution 


of the problem of two-level control is given by the formula 

|gCTl| S = IPPjCOl - 11 _1 IPjCTl + PjCOI - 11, C171 


where PjCOI and PjCTI are the input and the output values of the 
first level population, respectively. Since the evolution is 
unitary one has P^Ctl + PpCtl * I* 
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4. POPULATION CONTROL OF A RESONANTLY DRIVEN THREE-LEVEL SYSTEM 


Let us consider the coherent dynamics of a three-state system 
with an arbitrary level configuration simultaneously driven in 
resonance by three lasers with arbitrary and different time 
dependences of their amplitudes. Choosing the phases in an 
appropriate way we can write the system Hamiltonian in the matrix 
form with the following nonvanishing elements 


H 12 * H 21 " h 12 CU » H 23 * H 32 = h 23 m * H 13 “ " H 31 * -ihjgCtl. 

If we write this matrix in terms of the generators of the standard 

representation of the SUC21 group, then the time-evolution matrix 
is simply the three-dimensional representation of this group 171 


U 


m 


' u 

-V?uv* 

Cv*} 2 


V2uv 

|u| 2 -|v| B 

-V2u*v* 


v 2 


V2lI*V j 


Cl 81 


£ 

where u e g and v = - C2ig - hQgl/2h. The governing equation for 
the driven three-level system is given by 

g - Ifh i2 -ih 23 5fh 12“ ih 23 ) ” 1 + ih 13^ g + + hp.3) 9 a 0. 

The complete formal solution of the problem of three-level 
population control is given by the formula 

CPjCOl-PjpCOllCP^Tl-PgCOll- CP^Ol-PgCOlXP^CTl-PpCO}} 

|gCT3| 2 - . 

CP 1 C01-P 3 C011CP 1 C01+P 3 C01-2P 2 C011 

The useful relation exists between two variables u and v 

| uC 1 3 | 4 - | vC 1 3 | 4 = [PjCOl - PgCOir 1 IPjCtl - P 3 Ctll. 

Using the eqs. 10 and C181 we can now easily find simple exact 
conditions for complete transfer of population from one level to 
another and for complete return of population to an initial level. 
Let at t*0 level 1 is fully populated, i.e.PjCO}* 1, P^COl* P 3 C01* 
0, then the population is transfered completely to the level 3 at 
time t, i.e. P 3 Ctl* 1, P|Ct1 » PpCtl* 0, if gftl * 0. For the same 
initial conditions the population returns to the first level at 
time t, i.e. PjCt) * 1, PpCt}= P 3 Ctl* 0, when |gftl | * 1. 
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5. POPULATION DYNAMICS OF A RESONANTLY DRIVEN FOUR-LEVEL SYSTEM 

The resonant Hamiltonian of a four-level system simultaneously 
driven by six lasers with arbitrary and different time dependences 
of amplitudes is given by the eq. 9 with N*4. If the phases are 
choosen in such a way that all the Rabi frequencies C83 are purely 
imagine ones then the four-level Hamiltonian generates the 
following 4x4-matrix realization of the S0C43 algebra 

32 13 21 23 31 12 14 24 34 41 42 43 

R, . r 2 - r 3 -Rj *-r 2 -r 3 • Jj - j 2 ■ J 3 — jj -1 2 -j 3 ■ 1, 

with zero remaining matrix elements and with the following 
commutation relations 

IR j> V * e ju R i> IJ j- J k> * e jki R i> ,R j> V * e jki J r 

One can rearrange this basis by such a way to obtain from the 

matrices R and J. two mutually commuting sets of matrices A + 

and B . J with the commutation relations C12). Now we can rewrite 

o, +, - 

our Hamiltonian in the form 

HCU = Z V 13 A » + V tJ B n. ’ C193 

m=o, +, - 

where the following short-hand notations are introduced 

a o = iCh 12 " h 34 } ' 2a - E ih 23 + h 13 “ lh 14 + h 24» 

2a + e ih 23 - h J3 -ih J4 - hg 4> b Q e iChj 2 + h 34 ), C201 

2b_ = ih 23 + h 13 +ih J4 - h 24 , b + = ih g3 - h 13 +ih 14 + h 24 . 

Since HCt3 exactly equals to the sum of the two SUC2) 
Hamiltonians, commuting with each other, we may use all the results 
obtained for the SUC21 dynamical symmetry. Thus the system 
posesses SUC2I ® SUC21 dynamical symmetry [81 and the evolution 
operator can be written in the factorized form U * U g U f , where 

Ug * exp CgQ AqD exp Cg_ A_3 exp Cg + A + 3, 

U f * exp Cf Q B q 3 exp Cf_ BJ exp Cf + B + 3 C213 
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The governing equations for our four-level system take the form 

• • • 1 * * O 

g - a + a” g + Cl/2)Ci ay - iaQ a + a” + 2a.a + + Cl/2)ay)g ■ 0, 

f’ - b + b~* f + Cl/2Xib 0 - ib Q t^b” 1 + 2b_b + + Cl/2)by)f ■ 0, 

with the variables g = exp Cg 0 /2) and f = exp Cf Q /2) respectively. 
In an explicit form the evolution matrices are the following 

g Reg g Img Reg_ -Img_ 

-g Img g Reg Img. Reg. 

-Reg_ -Tmg_ gReg -g Img 

Img. -Reg. g Img g Reg 

f Ref f Imf Ref. 

-f Im f f Ref Imf. 

y a j.-l -Ref. -Imf. f Ref 

-Imf. Ref. -f Imf 

where g. ■ g C2ig - ay g)/2a + and f. * f C2if - by f)/2b + . 

Let at t * 0 the first level is occupied with probability one, 
i.e. PjCOl* 1, PgCOl* P^CO)* P^CO) * 0. Then the exact condition 
for complete transfer of population to the upper level is given by 

|g _1 CRef Img. + Reg. Imf) - f _1 CReg Imf + Ref. Img.) | 2 ■ 1 C22) 

and the condition for the first level to be completely depleted of 
population is 

Reg Imf - Img Imf - Cgf)”* CReg. Ref. - Img. Imf.) * 0. C23) 
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